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Abstract. 

We investigate the system size dependence of jet- quenching by analyzing transverse 
J^ • momentum spectra of neutral pions in Au+Au and Cu+Cu collisions at yS^N =200 



GeV for different centralities. The fast partons are assumed to lose energy by radiating 
fSJ ' gluons as they traverse the plasma and undergo multiple collisions. The energy loss 

per collision, e, is taken as proportional E (where E is the energy of the parton), 
proportional to \/E, or a constant depending on whether the formation time of the 
gluon is less than the mean path, greater than the mean free path but less than the 
Q . path length, or greater than the path length of the partons, respectively. NLO pQCD is 

^ I used to evaluate pion production by modifying the fragmentation function to account 

for the energy loss. We reproduce the nuclear modification factor Rjs^js^ by treating e as 

the only free parameter, depending on the centrality and the mechanism of energy loss. 

^^ ' These values are seen to explain the nuclear modification of prompt photons, caused 

by the energy lost by final state quarks before they fragment into photons. These also 
l/^ \ reproduce the azimuthal asymmetry of transverse momentum distribution for pions 

^O ' within a factor of two and for prompt photons in a fair agreement with experimental 

: . data. 
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1. Introduction 

Exploration of tlie properties of quark gluon plasma (QGP) produced in relativistic 
heavy ion collisions represents a major theme of modern nuclear physics research. The 
formation of QGP is heralded by numerous signatures, e.g., jet-quenching [HE], elliptic 
flow of hadrons [31 Hj, quark number scaling of elliptic flow of hadrons[5], radiation 
of thermal photons [SI El El E], and suppressed production of J/V^ [101 E], etc. The 
celebrated phenomenon of jet-quenching has its origin in energy loss suffered by high 
energy quarks and gluons as they traverse the QGP, colliding with other partons and 
radiating gluons. It leads to a suppressed production of hadrons having large transverse 
momenta as compared to the case for NN collisions at the corresponding centre of mass 
energy per nucleon. It is measured in terms of the nuclear modification factor Raa^ 
given by: 

T} ( V, _ d'^Np^p^{h)/dpTdy .. 

where b is the impact parameter and TAA{b) is the nuclear overlap function for impact 
parameter b. 

The jet-quenching has several other manifestations. Thus, for example, in a non- 
central collision the jets moving in and out of the reaction plane would cover differing 
distances inside the plasma, lose differing amounts of energy, and would be quenched to 
a different extent. This will lead to an azimuthal anisotropy of momentum distribution 
of hadrons which does not have its origin in the so-called elliptic flow [12]. 

Next consider a high energy quark which is produced in a hard scattering. A 
photon fragmented off this quark contributes to the prompt photon production at high 
transverse momenta. The quark may lose energy before its fragmentation and thus 
this photon production may get suppressed [131 HI] • Following the same arguments as 
above, this may also lead to an azimuthal anisotropy in photon production for non- 
central collisions [15]. The jet-photon conversion [7] and jet-induced bremsstrahlung [8] 
lead to very small contributions with opposite signs for the azimuthal anisotropy of the 
photon productions [T5] and are not considered in the present work. 

A new dimension to the azimuthal asymmetry and the jet-quenching has been 
added by the realization that the event by event fluctuations of the initial conditions 
may affect the elliptic flow of the flnal state hadrons p^. It is not very clear these 
seriously affect the azimuthal variation of the transverse momenta, which arises from 
jet-quenching [17] as well. 

It is expected that the magnitude of all these effects will depend on the dynamics of 
the collision and the properties of the QGP. This makes jet-quenching a powerful probe 
of the QGP. 

While there are quite a few detailed studies on the jet-quenching at RHIC and 
LHC energies ^8], in the present work we continue the study reported in Ref. [T3l [22] , 
where the Wang, Huang, and Sarcevic model [23] was used to study the evolution of 
the mechanism of energy loss at RHIC and LHC energies in central collision of heavy 
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nuclei. The model treats the medium produced in the collision as a plasma of length L 
through which the partons move, undergoing multiple collisions, losing energy through 
radiated gluons, before fragmenting to form hadrons. The only parameters which then 
enter the calculations are the path-length- L, the mean-free path- A, and the energy loss 
per collision- e. 

The energy loss per collision is then assumed to take a form as proportional to the 
energy of the parton- E^ \fE^ or a constant; depending upon the formation time of the 
radiated gluons, which is less than the mean free path, larger than the mean free path 
but less than the path length L, and greater than the path length L, respectively. The 
three regimes are called the Bethe-Heitler regime, the LPM regime, and the complete 
coherence regime (as the entire medium radiates as a whole) [2^. The formation time 
(or length) of the radiated gluon is ~ w/fcf^, where u is the energy of the gluon and kx 
is its transverse momentum. 

In the present work we explore the limits of the above approach by performing a 
detailed and systematic study of medium modification of production of neutral pions and 
prompt photons having large transverse momenta as a function of centrality and system 
size, by analyzing the corresponding data at 200 GeV/nucleon for collisions of gold nuclei 
and copper nuclei. We first seek an accurate description of the medium modification 
parameter Raa for different centralities and pt-, using the three approaches for the 
energy loss mechanisms, by varying the coefficients of energy loss, while accounting 
for the variation of the path length L. We find that the hadrons having low pt (< 5 
GeV/c) fall into the Bethe-Heitler regime, those having intermediate pt (5 GeV/c < pt 
< 10 GeV/c) fall into the LPM regime, while those having large pt > 8-10 GeV/c are 
covered by the complete coherence regime. We find a systematic decrease of the energy 
loss per collision as we go to more peripheral collisions from central collisions, as one 
might expect. 

We then use the same parameters to determine the azimuthal anisotropy coefficient 
V2{pt) for neutral pions having large pt- We find that in general, the experimentally 
obtained V2{pt) is lower than those predicted by our calculations, though the trend 
is generally reproduced. Part of this discrepancy may have its origin in the uniform 
densities for the nuclei, used in the present work. 

Next we calculate the prompt photon production as well as vlikx)- Again we get 
reasonable description of the experimental data. Further improvements can perhaps 
be made by relaxing the condition of equality of mean-free path and energy loss per 
collision for quarks and gluons assumed while obtaining a description for R\ a . This 
would, however, increase the number of parameters. 

We give our formulation for particle production, multiple scattering, and energy 
loss as applied to pion production at 0{a1) in Sec. 2. The centrality dependence of 
nuclear modification of neutral pions and parton energy loss (dE/dx) for Au-Au and 
Cu-Cu collisions at the highest RHIC energy are discussed in Sec. 3. The results for the 
azimuthal asymmetry of high px neutral pions are shown in Sec. 4. In Sec. 5, we give the 
results for the nuclear suppression and azimuthal anisotropy of prompt photons which 
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arise due to the modification of tlie fragmentation contribution. Finally we summarize 
our findings in Sec. 6. 

2. Formulation 

2.1. Particle production in pp collisions 

In perturbative QCD, the inclusive cross-section for the production of a particle in 
proton-proton collisions is given by [251 ES] 

d^PTdy 



da^^^^ ^ /■ , /■ , r dz 

K/A[Xa,Q^) 

dcTab^c 



J2 dXa j dXb —Fa/A{Xa,Q'^)Fb/B{Xb,Q'^) 

a,b,c 



where Fa/Ai^jQ"^) is the parton distribution function (PDF) for the parton a and 
Fb/B{x,Q'^) is the PDF for the parton b, for the nucleons A and B respectively. Dc/c 
gives the fragmentation probability of parton c into a particle C. In case of a hadron; 
Dc/h is the fragmentation function evaluated aX z = Ph/Pc-, where z is the fraction of the 
parton's momentum carried by the hadron. In case of photon production, -Dc/7 gives 
the fragmentation probability of a photon fragmented off a quark with the momentum 
fraction z = p-y/pc- In addition, we have an extra term where photon is directly produced 
in the hard collision (0 = 7). In that case the fragmentation function reduces to 6{l — z). 
aab^c is the parton-parton cross-section, calculated for the leading order processes C(a^) 
such as: 

g + g -> g + g, 
g + g -> g + g, 

q + g-^ q + 9, 

9 + 9^9 + 9, 

(3) 



At the next-to- leading order, C(a^) we include subprocesses like: 

g + g-> g + g + fi- , 
g + g-^g + g + 5f, 
g + g' -> g + g' + fi-, 
g + g -^ g' + g' + 5f, 

9+9^9+9+9 



(4) 



The running coupling constant as(/i^), is calculated at next-to-leading order 



Un ( 6(153 -19A^/) In ln(/iVA2 



(33 - 2Nf) ln(/i2/A2) \^ (33 - 2NfY ln(/i2/A2) 



(5) 
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where /i is the renormahzation scale, Nj is the number of flavours, and A is the Aqcd 
scale. We have used the CTEQ4M structure function [27] and the BKK fragmentation 
function ^E\- We set the factorization, renormahzation, and fragmentation scales as 
equal to pt, though we have checked scale dependence on particle production in pp 
collisions using Q = 0.5pT and Q = 2.0pr as well. 
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Figure 1. (Colour on-line) A comparison of neutral pion yield in p+p collision 
measured by the PHENIX collaboration [H] at ^s^N = ^00 GeV with NLO pQCD 
calculations. 



We show the results of our calculations for the inclusive neutral pion yield for 
proton-proton collisions at the centre of mass energy of 200 GeV along with the data from 
the PHENIX measurements [29] in Fig. [1] A very good agreement with the experimental 
data is seen, as in similar calculations by other authors (see e.g., Ref. [3D]). The 
quantitative description of the data from pp collisions provides us with a reliable base- 
line to discuss nuclear modification of hadron production in nucleus-nucleus collisions. 



2.2. Model for parton energy loss 

In order to calculate the inclusive pion spectra from nucleus-nucleus collisions we have 
used the effects of shadowing and energy loss of partons. The nuclear shadowing 
accounts for the modification of the parton distribution function for nuclei compared to 
that for a free nucleon. The shadowing function is defined as: 

Fa/A{X,Q^) 



Sa/A{X,Q ) 



AF,/^(x,g2 



(6) 



where Fa/A{x,Q^) is the parton distribution for the nucleus A and Fa/N{x,Q'^) is that 
for a nucleon. We have used the EKS98 parametrization of nuclear shadowing obtained 
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by Eskola, Kolhinen, and Salgado [3T]. 

We use the simple phenomenological model, first used by Jeon et al. [13] to obtain 
a description of preliminary Raa at RHIC energies, to estimate parton energy loss. The 
key features of this treatment are: average energy loss per collision (e^) by the parton, 
the mean free path of the parton (Aa), and the average path length of the parton in the 
medium (L). The probability for a parton to scatter n times before leaving the medium 
can be written as: 

The multiple scattering and energy loss suffered by the partons then necessitates a 
modification of the fragmentation function Dc/h{z^Q^)- Following the model of Wang, 
Huang, and Sarcevic [23], the modified fragmentation function can be written as: 

N 



'^N n=0 



^n-^c/Zil^niQ ) + 2^ ^m^g/hi^myQ 



(8) 



^m ■'-^g/hV'^mi "^ 
m=l 

where zE^ = z'^E^ = z^{E^ — J27=o^a)^ ^m = ^^t/^T- ^ is the maximum number 
of collisions for which z'^ < 1, D^^,^ is the hadronic fragmentation function which gives 
the probability that quark or gluon would fragment into a hadron (in our case vr") and 
C% = Y^n=o Pain, L) . The first term gives the fragmentation probability of the leading 
parton a with a reduced energy (E^ — X^ILo ^D ^^^ the second term comes from the 
emitted gluon having energy e™. 

While one may easily ignore the azimuthal (0) dependence of the average path 
length L[b) for impact parameter b for very central collisions, it needs to be considered for 
non-central collisions. We use a simple approach, based on Glauber model, to evaluate 
the dependence of the average path-length on the azimuthal angle with respect to the 
reaction plane. Assuming uniform densities for the colliding nuclei, the average path- 
length for an impact parameter b and azimuthal angle can be written as: 

}} TAB{x,y;b)dxdy 

where x and y are the transverse co-ordinates for the point where the partons scatter 
to produce the jet(s) which, travelling at an angle with respect to the reaction plane, 
traverses the path length i{x, y, 0, b). Tab^x, y; b) = t^(x -|- 6/2, y)tB{x — 6/2, y) is the 
nuclear overlap function and Ia and t^ are the transverse density profiles of the two 
nuclei. An average of L(0; b) over (varying from zero to 27r) gives the average path 
length L{b) (see Fig.Ej). 

We closely follow the excellent review on jet-quenching by Baier et al [21], and 
assume that the collisional energy loss is quite small for light quarks and gluons 
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Figure 2. (Colour on-line) (Left panel) Azimuthal variation of the average path length 
traversed by a parton in collision of gold nuclei. The impact parameter for the upper 
curve is an average for 0-10% most central collisions and the lower one is for 40-50% 
centrality. (Right panel) The average path length vs. impact parameter for Au-|-Au 
system. 



compared to the energy energy loss due to radiation of gluons. Now one can define 
the formation time of the radiated gluon as: 

form ~ ITi ^ ) 

rhrp 

where u is the energy of the radiated gluon and kx is the transverse momentum. We 
would normally have a; S> fc^ where fcy ~ /^ is a typical momentum transfer in a partonic 
collision. The coherence length, /^qJ^, can now be written as: 



(^T)coh ^coh(^T) 
where N ^ is the number of coherent scattering centers. One can then write, 



N. 



I 



coh 



coh 



u 



u 



Ej 



11 



(12) 



where -^lpm = Aq/i^ is the energy parameter introduced to separate the incoherent and 
coherent radiation of gluons. 

In the Bethe-Heitler (BH) regime, the energy of the radiated gluon is less than 
-^LPM' ^^^ coherence length l^^^i ^ Aq, leading to an incoherent radiation from L/A^ 
scattering centers. Thus energy loss per unit length in the soft u limit becomes: 

(13) 

where N^ = 3 and E is the energy of the parton. We shall write £„ ~ kE for this case 

and determine k from the measurement of Rj^j^- 

For u > -E'lpm, we have a regime of coherent radiation, where the coherence length 

^coh ^^ greater than A^ but less than average path length L. Thus, ^^oh^^ ^^ centres 

of scattering radiate coherently leading to 

dE a.Nr r^ 

^LPM^ . (14) 



dE 




(Tn 




1 




dx 
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TT 


^c 
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dx 



vr A„ 
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We re-write this as Sa ~ yoE and determine the values of a from the measurement of 

Finally, when the coherence length is larger than L {Iqq\^ > L) the entire medium 
acts as one coherent source and energy loss per unit length becomes: 

- ^ = ^AT.Ml . (15) 

dx TV Xa 

This is denoted as constant energy loss {Sa = k) regime. We use the NLO code of 
Aurenche et al [25] adopted earlier by Jeon et al. [13] to account for the energy loss 
of the partons, by modifying the fragmentation function, for the calculations reported 
here. 

A more detailed derivation of the energy loss for this case leads to: 

_ ^ . ^nJJ^, . (16) 

dx A Xa ^ ^ 

where "v is Fourier transform of the normalized differential cross-section for parton-parton 
collisions for the appropriate momentum transfer scale (see Ref. [24j). The momentum 
transport coefficient q is then defined as: 

g = Y^i), (17) 

SO that we can write: 

which can be used to deduce the average momentum transport coefficient for any given 
centrality, for a comparison with other estimates for this. 

3. Results 

3.1. Centrality dependence of Raa for Au-Au collisions at RHIC 

We have calculated the nuclear modification factor Raa for neutral pions for Au-Au 
collisions at the top RHIC energy for six centralities using Bethe-Heitler (BH), LPM, 
and constant energy loss mechanisms, discussed above. We have kept the mean free 
path of gluons as well as quarks fixed at 1 fm and also assumed that their energy loss 
per collision is identical. Thus once we have accounted for the variation of the average 
path length L{b) with centrality, we are left with the energy loss per collision as the 
only adjustable parameter. 

The results for BH mechanism are shown in Fig. Ofor near-central (0-10%, 10-20%, 
20-30%) and mid-central (30-40%, 40-50%, 50-60%) collisions. We see a good description 
of the data for pt up to about 6 GeV/c. The necessary energy loss per collision is seen 
to slowly drop from 10% to about 8% of the energy of the partons as we move from near 
central to mid-central collisions. We also note that the slope of Raa changes around pt 
equal to 5 GeV/c, which suggests a possible change in the mechanism for energy loss 
for partons contributing to higher momenta [22] . 
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Figure 3. (Colour on-line) Nuclear modification of neutral pion production for Au+Au 
collisions at a/snn=200 GeV, using BH mechanism (see text). The experimental data 
are taken from Ref. i32i. 



This expectation is confirmed by results shown in Fig. HI where the so-called LPM 
mechanism of energy loss is seen to provide an accurate description of the nuclear 
modification function in the pt range of 6-10 GeV/c. Once again we note that the 
energy loss coefficient a decreases systematically as the collisions become less central. 
Thus a parton having an energy of 10 GeV is likely to lose about 1 GeV in the first 
collision in the most central collision considered here and about 0.7 GeV in the collisions 
having a centrality of about 50-60%. 

Finally, we see (Fig. [5]) that the medium modification function for hadrons having 



System size dependence.... 



10 



0-10%; Au+Au@200A GeV 



LPM; E=V(aE) 



-- a=0.110 
— a=0.120 
- a=0.130 



i-i. 




-- a=0.100 

10-20%; Au+Au @ 200 AGeV — a=0.110 

- a=0.120 

LPM; E=V(aE) 








I ' ' ' I ' ' ' I 



-- a=0.090 

20-30%; Au+Au@200A GeV — a=0.100 

- a=0.110 



LPM; E=V(aE) 







2 4 6 8 10 12 14 16 

p^ (GeV/c) 



2 4 6 8 10 12 14 16 

p^ (GeV/c) 



2 4 6 8 10 12 14 16 

p^ (GeV/c) 



- a=0.100 
— a=0.090 
-- a=0.080 



i?*^ 



==^3*T 



30-40%; Au+Au@200A GeV 
LPM; E=V(aE) 



I ' ' ' I ' ' ' I 



C f*Tf 



iMfe 



40-50%; Au+Au@200A GeV 
LPM; E=-V(aE) 



-- a=0.060 
— a=0.070 
- a=0.080 








-- a=0.040 

— a=0.050 

50-60%; Au+Au@200A GeV -0=0.060 



LPM; E=V(aE) 



2 4 6 8 10 12 14 16 

p^ (GeV/c) 



2 4 6 8 10 12 14 16 

p^r (GeV/c) 



2 4 6 8 10 12 14 16 

Pt (GeV/c) 



Figure 4. (Colour on-line) Nuclear modification of neutral pion production for Au+Au 
collisions at y/sj^=200 GeV using LPM mechanism. The experimental data are taken 
fromRef. 32;. 



Pt > 8 GeV/c is best described by the mechanism assuming a constant energy loss per 
colhsion. The energy loss per collision is seen to decrease from about 1.4 GeV for the 
most central collisions to about 1 GeV for the collisions having a centrality of 50-60%. 

The transition from the BH regime to the LPM seen above and noted earlier [22] 
at pt ~ 5 GeV/c should not come as a surprise to us, since for the mean free-path 
A ~ 1 f m and (k^) ~ 1 (GeV/c)^ per collision, -E'lpm is about 5 GeV, which roughly 
separates the BH and LPM energy loss regimes for all centralities 

The results for these three figures are summarized in Fig. El Several interesting 
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facts emerge. 

Let us first look at the results for the BH regime operating for p^ up to about 5-6 
GeV/c. We note as we go to more peripheral collisions the data, even for lower pr, 
are in a better agreement with our calculations. This, we feel, is due to the reducing 
importance of radial flow (which gives a transverse kick to the partons/hadrons) in 
peripheral collisions. 

The LPM regime operating over the pt window of 5-10 GeV/c and the complete 
coherence regime operating at higher pt are seen to correctly describe the details of the 
changing curvature of the nuclear modification factor. 

We add that a look at Fig. |l]may prompt one to conclude that the description using 
the LPM mechanism is also working reasonably well till the largest pt considered here. 
However a closer look at the Figs. H] and |5] reveal that the LPM description mostly misses 
the data for larger pt values while the description using the constant energy loss per 
collision correctly reflects the curvature of the data till the largest pt- A more accurate 
data extending up to even larger px would be very valuable to settle this question more 
firmly. 

Even though we have identified the regions of the applicability of the mechanisms 
of energy loss, purely on the basis of a good description of the experimental values of 
the nuclear modification factor, it is heartening to note that the separation of the BH 
and the LPM regimes is very close to 5 GeV/c expected by us, as discussed above. We 
recall that the LPM and the complete coherence regimes differ by the formation time of 
the gluons, which is larger than the mean free path but smaller than the path-length for 
the former, while it is larger than the path-length (and of-course the mean free path) for 
the later. This is reflected in a slight change in the curvature of the nuclear modification 
factor, around px ~ 8-10 GeV/c. 

In brief, we have noted that the energy loss parameter in the Wang, Huang, and 
Sarcevic model for energy loss of partons in relativistic heavy ion collisions at RHIC 
can be tuned to obtain an accurate description of the nuclear modification factor Raa 
for different ranges of the transverse momenta. This observation is further strengthened 
by Fig. [7] where the variation of Raa integrated over px ranging from 6 to 9 GeV/c 
is shown as a function of number of participants and compared with the data from 
PHENIX collaboration [331. 



3.2. Centrality dependence of Raa for Cu-Cu collisions at RHIC 

In the next step we apply the above treatment to the collisions of copper nuclei at the 
same centre of mass energy of 200 GeV/nucleon. These collisions present an interesting 
system for the study of jet-quenching which complements as well as supplements the 
results from the collision of gold nuclei discussed earlier. It is well known that the 
more central collisions of copper nuclei have number of participants similar to those 
for mid-central collisions of gold nuclei. Thus a comparison would give us results for 
two systems, involving similar number of participants, one of which has nearly central 
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Figure 5. (Colour on-line) Nuclear modification of neutral pion production for 
Au+Au collisions at y/sj^=200 GeV assuming a constant energy loss/collision. The 
experimental data are taken from Ref. |32) . 



collisions and the other one has a large eccentricity. 

Proceeding as before, we give our results for the nuclear modification factor Raa for 
Cu+Cu collisions for centralities varying from 0-10% to 30-40%, in Figs. [Hl-fTOl using 
the BH, the LPM, and the constant energy loss mechanisms and compare with the 
measurements by PHENIX collaboration [34J. (See also Ref. [35]). 

From Fig. [HI we see that the energy loss mechanism suitable for the BH regime 
provides an accurate explanation of the data for all the centralities under consideration 
for p2- < 6 Q€\[ jc. 
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Figure 6. (Colour on-line) Nuclear modification of neutral pion production for Au+Au 
collisions at -y/SNN=200 GeV for the three energy loss schemes for different pr regimes. 



The energy loss mechanism denoted as LPM is seen to describe the nuclear 
modification function over the pt range of 6-10 GeV/c and even beyond for less central 
collisions Fig. |9l 

And finally, the results for the mechanism admitting a constant energy loss per 
collision is seen to work well for pions having px > about 6-8 GeV/c (see Fig. [TU]) . 

We note that the energy loss co-efficient for the most central collision of copper 
nuclei is close to that for the collision of gold nuclei for the 50-60% centrality, for all the 
three mechanisms. We also note the possible separation of the BH and LPM regimes 
around pt equal to 5-6 GeV/c, as before for gold collisions, though the demarcation of 
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Figure 7. (Colour on-line) Centrality dependence of Raa of neutral pions calculated 
for LPM mechanism of energy loss. The data are from the PHENIX collaboration [33] 
for Au+Au collisions at J's^^ — 200 GeV. 



the LPM and the constant energy loss mechanisms is not as sharp for the copper nuclei 
as for the gold nuclei, especially for the less central collisions. This may have its origin 
in smaller L for copper nuclei. 

A figure similar to Fig. [6] can also be prepared for the collisions involving copper 
nuclei. We have not given it here, for the want of space. 

3.3. Centrality dependence ofdE/dx 

The results for the largest values of transverse momenta are summarized in Fig. [TTl 
Taking the case of constant energy loss per collision, —dE/dx = e/X, we see that the 
rate of energy-loss necessary to explain the suppressed production of hadrons at large 
Pt varies as (L), the average path length of the partons in the medium, both for Cu+Cu 
and Au+Au collisions at the centre of mass energy of 200 GeV/nucleon. This empirical 
result confirms the conviction that for the range of transverse momenta covered, Pt > 
8-10 GeV/c, the concerned partons interact with the medium as a whole, and the energy 
lost by the partons, AE oc L^. This is different from the AdS/CFT description where 
the radiated parton cloud is brought on-shell by the drag of the strongly coupled plasma 
and AE oc L^ [36j, though some moderation of the L dependence is expected for an 
expanding plasma. 

We also see that even though the slopes of the results for Cu+Cu and Au+Au 
collisions are similar, the energy loss for a given (L) for Au+Au collisions is about 
40-60% larger than that for Cu+Cu collisions. 

We can use the Eq. [18] to obtain the value of the momentum transport coefficient, 
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Figure 8. (Colour on-line) Nuclear modification of neutral pion production for Cu+Cu 
collisions at a/S]\j]\j=200 GeV for BH mechanism of energy loss along with the results 
from the PHENIX collaboration El- 



g from the above. We find that the g varies from 0.25 GeV^/fm for the collisions having 
0-10% centrahty for Au+Au system at 200 AGeV to 0.39 GeVVfm for the 50-60% 
centrality. The decrease for more central collisions can be understood by noting that 
the systems produced in such cases would be at higher temperatures and thus much 
more dense than for peripheral collisions. This would then lead to a smaller suppression 
of the radiation due to the LPM effect. The faster rise with the path- length (L^) more 
than compensates this decrease, finally giving a large energy loss and larger suppression 
of high momentum particles for more central collisions. 

On the other hand, for the case of Cu+Cu system, we get q ^ 0.18 GeV^/fm for all 
the centralities considered by us. We note that it is about half of what we got for Au+Au 
system. The near identity of q for all the centralities for Cu+Cu system is related to a 
smaller variation in the path-length for them. We also note that the Au+Au system is 
more efficient in degrading the jet-energy. 

We add that various authors have reported widely differing values of q. Thus, 
Wiedemann and Salgado [IB] have reported a value of 5-10 GeV^/fm, while Gyulassy, 
Levai, and Vitev [18] report a value in the range of 0.35-0.85 GeV^/fm. Arnold, Moore, 
and Yaffe [18], on the other hand, suggest a value of about 2 GeV^/fm, in contrast to a 
value of 3-4 GeV^/fm reported by Majumder [TSl [201 EI] • It is felt that a part of this 
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Figure 9. (Colour on-line) Nuclear modification of neutral pion production for Cu+Cu 
collisions at ^/sj^=200 GeV for LPM mechanism of energy loss and compared with 
PHENIX results |M1- 



difference may be due to different physical attributes of the evolving system over which 
the average is taken. 



4. Azimuthal anisotropy of high p^ neutral pions 

We have seen that the phenomenological model of parton energy loss, which accounts 
for only the induced gluon radiation of the fast partons, successfully describes the 
suppression of neutral pion spectra at large px at different centralities for Au+Au and 
Cu+Cu collisions at J^^^ = 200 GeV. In case of non-central collisions, the path-length 
will have an azimuthal variation over the transverse plane (see Fig. |2]). In what follows, 
we study the azimuthal dependence of the suppression of production of hadrons in terms 
of the azimuthal anisotropy of their transverse momenta [12]. 

We re-emphasize that this anisotropy is different from the elliptic flow seen for 
hadrons at lower p^, which is explained using hydrodynamics which converts the spatial 
anisotropy of the initial state into the azimuthal anisotropy of the transverse momenta 
of the hadrons ^. The extent of validity of hydrodynamics for larger p^ is not very 
clear. The success of the recombination model [5] suggests that the fragmentation as a 
mode of hadronization may be valid for pt ^ 3-5 GeV/c. Thus the following results 
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Figure 10. (Colour on-line) Nuclear modification of neutral pion production for 
Cu+Cu collisions at ^/snn=200 GeV for constant energy loss mechanism compared 
with PHENIX results El. 



should be taken as indicative of anisotropy which can arise due to medium modification 
of fragmentation function due to energy loss of partons. 

The differential azimuthal anisotropy is measured in terms of the parameter V2{pt), 
which is the second Fourier coefficient of the azimuthal distribution of hadrons in the 
reaction plane: 

/o" d(j) cos{2(j))dN/d'^pTdy 



V2{P7 



Jo"" d<i) dN / d^prdy 



(19) 



We have calculated the azimuthal anisotropy coefficient V2{pt) of neutral pions 
for pt > 2 GeV/c for Au+Au colhsions at a/ ■5^^=200 GeV for the six centralities 
mentioned earlier. The energy loss per collision is taken from the earlier analysis and 
dependent distribution of the pions is calculated by incorporating the dependence of 
the path-length L. 

The results of our calculations are displayed in Fig. [12] along with the data 
available from PHENIX collaboration [33] . We show the theoretical results only over 
the pt window where the corresponding mechanism was found to describe the nuclear 
modification data, earlier (see Fig. [6]). 

Recalling that the BH mechanism was found to provide a good description of the 
nuclear suppression, up to pt equal to 5-6 GeV/c, while the mechanisms for the LPM 
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Figure 11. (Colour on-line) Variation of dE/dx with {L) for Cu+Cu and Au+Au 
collisions at 200 AGeV for the constant energy loss per collision regime obtained in 
the present work. 



and the constant energy loss per collision provided description for higher px-, the results 
are quite interesting. First of all, even for pr ^ 6 GeV/c while the BH mechanism leads 
to a V2{pt) increasing with px, the experimental results show a contrary behaviour. 
It could be that a hydrodynamic expansion of the system under viscosity affects the 
hadrons having px up to 5-6 GeV/c [37]. The other two mechanisms, though, give a 
behaviour similar to the experimental data, in the region of transverse momenta where 
they are applicable. The theoretical values are larger by up to a factor of two in the worst 
cases of more central collisions, though they agree reasonably well with the experimental 
values for the least central collision considered here. 

It is interesting to see the centrality dependence of the integrated azimuthal 
anisotropy coefficient in the px range of 6-9 GeV/c. For this purpose, we have used the 
differential azimuthal anisotropy coefficient V2{pt) for the LPM mechanism obtained 
earlier (see Fig JT3|l . We see that in general our calculations reproduce the trend of 
the variation of V2 with the number of participants, though the theoretical values are 
larger by about a factor of two (see Ref. [38], for similar results). One short-coming of 
our calculation immediately comes to mind- the description of the nuclei as having a 
uniform density. A Woods-Saxon density profile for the colliding nuclei would reduce 
the difference in the path-lengths for the partons travelling in the reaction plane and 
travelling in a plane perpendicular to it, for example, and thus reduce V2. We postpone 
this study to a future work where additionally, the evolution of the medium will be 
taken into account. 

We also give our prediction for V2{pt) for Cu+Cu collisions at 200A GeV for two 
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Figure 12. (Colour on-line) The differential azimuthal anisotropy coefRcient V2 of 
neutral pion calculated using BH, LPM, and constant energy loss mechanisms for 
Au+Au collisions at ■>/snn=200 GeV. The experimental data are from the PHENIX 
collaboration [33]. 



typical centralities (Fig. [T^ . A behaviour similar to that for Au+Au collisions is seen. 
5. Nuclear suppression and azimuthal asymmetry of prompt photons 



Prompt photons, as indicated earlier, originate from (a) quark-gluon Compton scattering 
(g + 5^ — > g + 7), (b) quark-anti-quark annihilation {q + q ^ g + 'j), and (c) collinear 
fragmentation of a high energy quarks {q — ?► g + 7) following a hard collision. Of 
these the third process is affected by the energy loss suffered by the quark before its 
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Figure 13. (Colour on-line) The centrality dependence of integrated V2 of neutral 
pions calculated for LPM mechanism of energy loss for Au+Au collisions J's^^ = 200 
GeV. The data points are taken from PHENIX collaboration [55] . 
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Figure 14. (Colour on-line) The differential azimuthal anisotropy coefficient V2{pt) 
of neutral pion calculated using BH, LPM, and constant energy loss mechanisms for 
Cu-|-Cu collisions at y^sj^=200 GeV. 



fragmentation [131 El- It will also lead to an azimuthal anisotropy of the momentum 
distribution of hard photons [13] . 

Once again, we use the NLO pQCD code of Aurenche et al [26], suitably modified 
to account for the energy loss suffered by the quarks before they fragment into photons. 
As a first step we show our results for proton-proton collisions at 200 GeV (Fig. [T5|l 
using the fragmentation, factorization, and the renormalization scales as equal, with 
Q = A;t/2, kx and 2kT. We have used the CTEQ^M parton distribution function and 
the BKK fragmentation function for these calculations as well. We see that our results 
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are consistent with the data from PHENIX [39] and also with earher calculations (see 
Ref. [32]) along similar lines. For the rest of the calculations we have used a common 
scale Q = kx- 
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Figure 15. (Colour on-line) A comparison of production of prompt photons in p+p 
collision measured by the PHENIX collaboration [39] at ys^N ^ '^^^ ^^^ '^^^'^ ^L*-* 
pQCD calculations. 



While performing the calculations for the nucleus-nucleus collisions we include 
nuclear shadowing and energy loss suffered by the quarks as before. We also properly 
account for the isospin of the participating nucleons, which is essential while performing 
calculations for photons. The nuclear modification factor R\ a is then obtained in a 
manner similar to Eq. [TJ 

We show our results for the nuclear modification of hard photon production using 
the energy loss parameters obtained earlier along with the experimental results for 
collision of gold nuclei at 200A GeV in Fig. [161 We give the theoretical curves over 
the regions where the nuclear modification for pion production is in agreement with the 
experimental data (see Fig. (6]). 

We see that our results are in fair agreement with the data for kj- beyond 10 GeV/c 
where the mechanism of prompt photon production included in the present work is 
expected to dominate [H]. We do realize that for lower kx other mechanisms like jet- 
conversion [7j, induced bremsstrahlung [S] and thermal production [Hj will contribute. 

Next we calculate the the differential azimuthal anisotropy coefficient vJikr) for 
hard photons for Au-|-Au collisions at 200 A GeV. Recall that for lower kx (< 5-6 
GeV/c), the azimuthal anisotropy of the direct photons carries valuable information 
about the development of the elliptic flow of the plasma and momentum anisotropy of 
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Figure 16. (Colour on-line) Nuclear suppression of hard photons calculated using 
BH, LPM, and constant energy loss per collision for Au+Au (0-10%) collisions at 
^S]y]\^=200 GeV. The data points are taken from Ref. [4^ 
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Figure 17. (Colour on-line) The differential azimuthal anisotropy coefhcient vl of 
direct photons calculated using LPM scheme of parton energy loss for Au-|- Au collisions 
at ysj^=200 GeV along with experimental results from the PHENIX [42] (left panel) 
and the STAR ,43] (right panel) collaborations. 



the thermalized partons 

In Fig. [T7] we have shown the results for v^ikx) calculated using LPM and the 
constant energy loss per collision mechanisms of parton energy loss for two centrality 
classes (20-40% and 10-40%) and compared them with the recent results from the 
PHENIX |l2] and the STAR [43j collaborations. We make no further adjustments 
of the energy loss parameters obtained earlier. 

We see that the azimuthal anisotropy of direct photons for kx > Q GeV/c is 
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reproduced reasonably well from our calculations. The results of Figs. [16] and [17] come 
a pleasant surprise as we recall that while calculating the energy loss of the partons 
and the consequent suppression of hadrons at large pt, we did not distinguish between 
quarks and gluons. The photon production, however, is sensitive only to the energy loss 
suffered by quarks. One possible reason for this could be the dominating contribution of 
quarks to hadrons as well photons having large transverse momenta, which is then rightly 
sampled by our procedure. One could use this procedure to first fix the energy loss per 
collision for quarks and then proceed to determine the energy loss per collision for gluons 
by fitting the hadronic spectra. This would, however, require a much more accurate and 
extensive large transverse momentum data for photons. A further availability of isolated 
photon spectra along with the (total) single photon spectra will help us to subtract 
out the Compton plus annihilation contributions and increase the sensitivity of these 
calculations to energy loss suffered by quarks. 

6. Summary and Discussion 

We have used a simple and transparent model of multiple scattering and energy loss 
for partons passing through quark gluon plasma produced in relativistic collisions of 
heavy nuclei, before they fragment into hadrons or photons. Taking the centrality 
dependence of the average path-length into account and taking energy loss per collision 
as an adjustable parameter, we have obtained a description of hadronic suppression in 
Au+Au and Cu+Cu collisions at 200A GeV. The energy loss parameters are found to 
vary systematically with the centrality and the dE/dx for pt > 6-8 GeV/c is found 
to vary linearly with the average path-length (L), in agreement with the predictions of 
Baier et al. [21]. 

The same parameters are then used without any further adjustments to obtain the 
azimuthal anisotropy for high p^ pions. We find it to be about a factor of two larger than 
the experimental values. Next, we have calculated the suppressed production of hard 
photons and their azimuthal anisotropy, again with out changing any of the parameters. 
These are in reasonable agreement with the experimental findings. We have mentioned 
that in principal, the calculations could be improved by incorporating diffused density 
distributions for the colliding nuclei, instead of the uniform densities used in the present 
work. In that case, it remains to be seen if the f 2 could be used to further constrain the 
energy loss parameters. 

The success of the approach, considering that the dynamics of evolution, the 
temperature and flavour dependence of some of the parameters and a more detailed 
accounting of various interferences are not considered here, once again confirms the 
inadequacy of -R^A ^^ seriously distinguishing various details of the dynamics of the 
plasma |21j . 

Some improvements are definitely in order. The use of uniform density for the 
colliding nuclei can be easily improved upon. This, as we argued, could improve the 
description of azimuthal anisotropy of hadron production at higher px- One may also 
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take differing values for the mean free patlis for quarks and gluons and the energy loss 
per collision for them. Note, however, that if we argue that dE/dx for gluons is twice 
that for quarks and also that their mean free path is half of that for quarks, then the 
energy loss for quarks and gluons per collision will be identical. 

Why does this simple approach work so well? We realize that the px distribution of 
jets produced in primary hard scatterings is a rapidly falling function of the transverse 
momentum. Partons having a transverse momentum p^ after undergoing multiple 
collisions and losing a momentum Ap^ will populate the momentum space occupied 
by partons having transverse momentum p^ — Apx while the later will shift to still 
lower momenta. Thus if we make a good estimate of Ap^, we are likely to get a good 
description of the nuclear modification. 

But is this approach really so simple? NLO pQCD and accurate structure and 
fragmentation functions, along with the applicability of pQCD at larger transverse 
momenta ensure an accurate base-line for p+p collisions. The influence of the dynamics 
of evolution is perhaps softened by the short life-time of the plasma at RHIC energies 
and the fact that most of the energy is lost in collisions at early times. 

Thus we feel that the present study, at the very minimum, gives a reliable average 
of the rate of energy loss of partons in quark gluon plasma. We also see that a 
simultaneous description of nuclear modification of photon and pion production in 
nuclear collisions, holds out the hope of getting the energy loss suffered by quarks and 
gluons independently. 

It will definitely be of interest to see if the centrahty dependence of dE/dx ( oc (L)) 
for partons having larger transverse momenta at LHC energies (2.76 ATeV and 5.5 
ATeV) continues to follow the trend seen for the data at RHIC energies for Cu+Cu and 
Au+Au collisions, seen in Fig. [TTl If confirmed at LHC energies as well, this will be a 
very valuable empirical affirmation of the suggestions of Baier et al [21]. This aspect is 
under investigation. 
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